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We further investigate center vortex percolation and Coulomb gauge remnant symmetry breaking in the SU(2) 
gauge-Higgs model. We show that string breaking is visible in Polyakov line correlators on the center projected 
lattice, that our usual numerical tests successfully relate P-vortices to center vortices, and that vortex removal 
removes the linear potential, as in the pure gauge theory. This data suggests that global center symmetry is not 
essential to the vortex confinement mechanism. But we also find that the line of vortex percolation-depercolation 
transitions, and the line of remnant symmetry breaking transitions, do not coincide in the SU(2)-Higgs phase 
diagram. This non-uniqueness of transition lines associated with non-local order parameters favors a straightfor- 
ward interpretation of the Fradkin-Shenker theorem, namely: there is no unambiguous distinction, in the SU(2) 
gauge-Higgs models, between a "confining" phase and a Higgs phase. 



I. INTRODUCTION 

Investigations of the confining force generally concentrate 
on gauge theories in which confinement is permanent; i.e. the 
linear potential increases without limit. Theories of this kind 
(with finite rank gauge groups) are all invariant under a global 
center symmetry, which can be expressed as 

Uo(x,t) — >zUo(x,t) all x, fixed t (1.1) 

in lattice formulation, where z 7^ 1 is an element of the (non- 
trivial) center of the gauge group. The unbroken realization of 
this symmetry is responsible for the vanishing of Polyakov 
line expectation values, and hence permanent confinement. 
For SU(N) gauge theories with this global symmetry, the po- 
tential between static color sources, in color group representa- 
tion r, depends only on the N-ality of representation r. While 
this fact is easily understood in terms of energetics/string- 
breaking arguments (e.g. a flux tube between adjoint sources 
can "snap" due to pair production of gluons), it also means 
that the string tension of a Wilson loop, evaluated in an en- 
semble of configurations generated from the pure Yang-Mills 
action (and therefore blind to the location of the Wilson loop), 
depends only on the N-ality of the loop representation. This 
leads to a rather profound conclusion: large-scale vacuum 
fluctuations — occuring in the absence of any external source 
— must somehow contrive to disorder only the center degrees 
of freedom of Wilson loop holonomies. The center vortex 
confinement mechanism (c.f. ref. 1 1| for a review) is the sim- 
plest proposal for how this type of disorder can occur. 

However, not all gauge theories of interest are invariant 
under a non-trivial global center symmetry Jl.lt . Examples 
include real QCD, and any other SU(N) gauge theory with 
matter fields in the fundamental representation of the gauge 
group. Another relevant example is G(2) pure gauge theory, 
whose center symmetry and first homotopy group are both 
trivial. 1 In these theories, the asymptotic string tension is zero, 



In SU (N)/Zn pure gauge theory, which has a trivial center and zero asymp- 
totic string tension, vortices and vortex fluctuations are no different from 



and at large scales the vacuum state is similar to the Higgs 
phase of gauge-Higgs theory. Such theories are examples of, 
rather than exceptions to, the general statement that confine- 
ment is dependent on the existence of a non-trivial global cen- 
ter symmetry. On the other hand, real QCD and G(2) pure 
gauge theory, as well as gauge-Higgs theory in some regions 
of the phase diagram, have a static quark potential which rises 
linearly for some interval of color source separation, and then 
becomes flat. We will refer to this situation as "temporary 
confinement" , reserving the term "permanent confinement" 
for theories which have a non-zero asymptotic string tension 
for color sources in the fundamental representation. 2 

In a theory with temporary confinement, the simple (and es- 
sentially kinematical) motivation for the center vortex mech- 
anism is lost. Then it is not obvious that the center vortex 
picture, which is motivated by the N-ality properties of the 
asymptotic string tension, is relevant. The relevance (or irrel- 
evance) of vortices to temporary confinement is a dynamical 
issue, which we would like to investigate via numerical simu- 
lation. 

The simplest case to consider is SU(2) gauge-Higgs theory, 
with the scalar field in the fundamental (j = 1 /2) representa- 
tion. Two of our previous articles dealt with this model. The 
first, written in collaboration with R. Bertie and M. Faber @|, 
showed that P-vortices percolate when the couplings lie in the 
temporary confinement region of the phase diagram, and cease 
to percolate in the Higgs region, where there is no linear po- 
tential at all. We also found that center-projected Polyakov 
lines, in the temporary confinement region, show evidence 
of color screening by the scalar field. This work did not, 



those of SU (N) gauge theory, but the relevant Zn symmetry is that of the 
first homotopy group 2]. 
2 In both cases, of course, the asymptotic particle states are color singlets. 
But this is also true in a Higgs phase, where the condensate screens any 
external charge. A similar effect occurs in electrodynamics, for electrically 
charged particles placed in a plasma or a superconductor (l||. One does 
not normally refer to electric plasmas and superconductors as confining 
systems; what is going on is charge screening. We believe it is useful to 
distinguish between this kind of screening of particle charge, and whatever 
physics lies behind flux tube formation and the linear static quark potential. 
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however, attempt to show that P-vortices in center-projected 
configurations actually correspond to center vortices in unpro- 
jected configurations, as none of our usual tests for that corre- 
spondence were employed. A second article, in collaboration 
with D. Zwanziger |4], considered the spontaneous breaking 
of a remnant global symmetry, which exists after Coulomb 
gauge fixing, in the SU(2) gauge-Higgs theory. A confining 
color Coulomb potential is associated with the unbroken real- 
ization of this remnant symmetry, and it was found that rem- 
nant symmetry was unbroken in the temporary confinement 
region, and spontaneously broken in the Higgs region. We 
did not check, however, whether remnant symmetry breaking 
and vortex depercolation occur at the same place in the phase 
diagram (although we assumed this to be true). The reason 
was that the position of the depercolation transition, found in 
ref. I2B, was determined for a gauge-Higgs theory with vari- 
able Higgs modulus, while the position of the remnant sym- 
metry breaking transition, found in ref. |4], was computed in 
the frozen modulus version of the theory. The present article 
is intended to fill in these gaps in our two previous articles. 



II. CENTER DOMINANCE 

We consider a gauge-Higgs theory with a frozen modulus 
Higgs field. For the SU(2) gauge group, the action can be 
written as 

S = p £ \7i[UUUW] +y£iTr[<j)t (*)<>(*+£)] 

plaq x.P 

(2.1) 

where (]) is SU(2) group-valued. This theory was first studied 
numerically by Lang et al. |5]; the phase diagram is sketched 
in Fig. ^ There is a line of first order transitions, but only 
one thermodynamic phase; any two points in the diagram can 
be connected by a path which avoids all non-analyticity in the 
free energy. The absence of a transition completely separating 
the diagram into a confinement phase and a Higgs phase was 
demonstrated analytically by Fradkin and Shenker, and Oster- 
weiler and Seiler, in refs. |6]. Nevertheless, below the transi- 
tion line lies a temporary confinement region, where the static 
quark potential rises linearly up to some screening distance, 
while above the line the theory is Higgs-like at all distances, 
and the static potential is nowhere linear. 

We would like to study center dominance inside the tem- 
porary confinement region, but close enough to the transition 
line so that the the screening effect of the scalar field is de- 
tectable numerically. For this purpose, we compute the expec- 
tation value of Polyakov lines at p = 2.2, on an L 3 x 4 lattice. 
At P = 2.2, the first order transition occurs at about y = 0.84. 
The quantity we measure is 
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FIG. 1: Schematic phase diagram of the SU(2) gauge-Higgs system. 
The solid line is a line of first-order phase transitions. 



at y = and on an l) x Lj lattice, we must find 



(23) 



while for explicitly broken center symmetry (y ^ 0) it must be 
that (P) has a non-zero limit at large volume. 

Our data for Polyakov lines on the unprojected lattice, at 
P = 2.2 and y = 0, 0.71, is shown if Fig. [2 for lattice sizes up 
to 20 3 x 4. The straight line is a best fit through the y = data, 
and errorbars for some data points are smaller than the symbol 
size. It is clear that the y = data is consistent with eq. (12.31 . 
and (P) extrapolates to zero in the infinite volume limit. At 
y = 0.7 1 the system is still below the first-order transition line, 
and in the temporary confinement region. It appears from the 
data that at this coupling, (P) has stabilized (atL = 14, 16,20) 
to a non-zero value of (P) « 0.034(1). So at p = 2.2, y = 
0.71, color screening of Polyakov lines by the matter field is 
detectable. 

The data at these same couplings, for Polyakov lines on the 
center projected lattice is displayed in Fig. (|3}. The center 
projected data tells exactly the same story as the data on the 
unprojected lattice: at y = 0, Polyakov lines tend to zero at 
large volumes, while at y = 0.7 1 screening is detected, and 
the Polyakov lines stabilize at (P) w 0.120(4). This aspect of 
center dominance in the SU(2) gauge-Higgs model was previ- 
ously found in ref. @|, for the variable modulus version of the 
theory. 

We can go on to calculate the correlator of center-projected 
Polyakov lines (P(x)P(x+R)) atp = 2.2, y = 0.71,ona20 3 x 
4 lattice. The data is shown in Fig.|4] The dashed line is a best 
fit to the data, for R > 2, by the function 



f(R)= c + ciexp[-4a/?] 



(2.4) 



where P(x) denotes the Polyakov line passing through the 
point {x,f = 0}. In the case of unbroken center symmetry, 



From the fit we findco = 0.0182, (7 = 0.211. Not surprisingly, 
co is quite close to the square of the VEV of the Polyakov line 
in center projection, which is (P cp ) = 0.12 on the 20 3 x 4 lat- 
tice. In this way we see string-breaking, due to the dynamical 
matter field, from Polyakov line data on the center-projected 
lattice. 
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FIG. 2: Polyakov line values, on the unprojected lattice at P 
2.2, y = and y = 0.71, on L? x 4 lattice volumes with L 
10, 12, 14, 16, 20 The straight line is a best fit to the y = data. 



FIG. 4: Polyakov line correlator (P(0)P(R)) on the center-projected 
lattice. 
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FIG. 3: Same as Fig.|2| for Polyakov lines on the center-projected 
lattice. 



consistent with the assumed correspondence of P-vortices and 
center vortices. 
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Although the data displayed in the previous graphs makes 
a good case for center dominance in Polyakov lines in gauge- 
Higgs theory (which is not a new result), there is still the ques- 
tion of whether P- vortex excitations on the center-projected 
lattice correlate with gauge-invariant observables on the un- 
projected lattice. At y > global center symmetry is broken, 
and the 't Hooft loop operator B(C) |7] which creates a thin 
center vortex would not only raise the action at the loop lo- 
cation, but also on some surface bounded by the vortex loop. 
The position of this "Dirac surface" is no longer a gauge ar- 
tifact. We can still identify P-vortices via maximal center 
gauge fixing and projection, but the correspondence of these 
P-vortices to center vortices on the unprojected lattice cannot 
be taken for granted. Our standard test for this correspon- 
dence is to see if W\ (C) /Wb(C) — > — 1 in the large-loop limit. 
Here W„(C) represents a Wilson loop, computed from unpro- 
jected link variables, with the restriction that the minimal area 
of loop C is pierced by n P-vortices on the projected lattice. 
The result of this test, for spacelike loops on a 20 3 x 4 lat- 
tice at P = 2.2, y = 0.71 is shown in Fig. |5] It is much like 
the result found for pure gauge theories, and seems perfectly 



FIG. 5: Ratio of "vortex-limited" Wilson loops. W\ (C) is evaluated 
for loops pierced by a single P- vortex; Wq(C) is evaluated for loops 
which are not pierced by any P-vortices. 



Finally, we compare the Creutz ratios of projected and un- 
projected spacelike Wilson loops at p = 2.2, y = 0.71, on 
the 20 3 x 4 lattice, and look for the effect of vortex removal. 
The relevant data is displayed in Fig. [5] We see that the pro- 
jected Creutz ratios are constant for any R > 1, as in the pure- 
gauge theory, at roughly %(R,R) tzs 0.21, and the values for 
the Creutz ratios on the unprojected lattice also appear to con- 
verge towards this value. Note that this value for the asymp- 
totic string tension is consistent with the value obtained from 
Polyakov line correlators on the projected lattice. The ef- 
fect of vortex removal is also shown in Fig. |6] Vortices are 
removed via the de Forcrand-D'Elia prescription [8], which 
consists of fixing to maximal center gauge, and multiplying 
each link variable by its center-projected value. This is the 
vortex-removed ensemble. We see that the Creutz ratios in 
this ensemble go to zero asymptotically, just as in the pure 
gauge theory. 
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FIG. 6: Creutz ratios in the gauge-Higgs theory for unprojected, 
center-projected, and vortex-removed lattices. 



III. SYMMETRY BREAKING AND VORTEX 
PERCOLATION 

The Fradkin-Shenker theorem |6] assures us that there is no 
phase transition which completely isolates the temporary con- 
finement region from a Higgs phase; at least, no such transi- 
tion could be detected by any local order parameter. But what 
about non-local order parameters? Perhaps thermodynamics 
is not the ultimate arbiter, and there really exists some qualita- 
tive difference between the temporary confinement and Higgs 
phases, characterized by symmetry-breaking, or by conden- 
sation of solitonic objects, which is only detectable via non- 
local observables. A relevant example is the Ising model in 
the presence of a small external magnetic field h. In that case 
the global Zi symmetry of the zero-field model is explicitly 
broken, and there is no thermodynamic transition from an or- 
dered to a disordered state. On the other hand, there is a sharp 
depercolation transition in the h > as well as the h = case; 
the line of such transitions in the temperature-/; phase diagram 
is known as a Kertesz line \9]. In the gauge-Higgs model the 
coupling y > breaks the global Zj_ symmetry, and it is pos- 
sible that a sharp vortex depercolation transition could serve 
to distinguish the temporary confinement and Higgs phases of 
the theory. 

An alternative proposal for distinguishing these phases is 
associated with symmetry breaking. We know from the 
Elitzur theorem that a local gauge symmetry can never be 
spontaneously broken. On the other hand, after Coulomb or 
Landau gauge fixing there still exists some global remnant 
of the local symmetry, and these global symmetries can be 
spontaneously broken. So perhaps center vortex depercola- 
tion and remnant symmetry breaking define a unique Kertesz 
line, which unambiguously separates the temporary confine- 
ment and Higgs phases of the gauge-Higgs model | lOjfl . An- 
other candidate symmetry for distinguishing the two phases, 
advocated by the Pisa group, is a certain dual (abelian) mag- 
netic symmetry 1 11]. This approach, in the non-abelian the- 
ory, also requires fixing to some gauge. 

Since both the identification of vortices, and the definition 



of remnant (as well as magnetic) symmetries entails the choice 
of a gauge, the associated order parameters are non-local 
(if expressed as gauge-invariant operators), and the Fradkin- 
Shenker theorem does not rule out non-analytic behavior in 
such observables. On the other hand, if the transition lines as- 
sociated with each order parameter do not coincide, then the 
claim that any of these parameters can be used to "define" 
confinement, in the absence of a non-vanishing asymptotic 
string tension, becomes less compelling. In this section we 
will report on our results for the Kertesz lines corresponding 
to vortex depercolation, and to remnant symmetry breaking in 
Coulomb gauge. 



A. Symmetry Breaking 

We begin by reviewing some points made in ref. |4]. First 
of all, there is a remnant symmetry in "minimal" Coulomb 
gauge, defined as the gauge with minimizes, on the lattice, 



/? = -££ReTr[C/*(*) 

x k=\ 



(3.1) 



Fixing to this gauge still allows the following "remnant" 
gauge transformations which are global in space, but local in 
time: 

U k (x,t)=g(t)U k (x,t)g\t) , U (x,t)=g(t)U (x,t)g\t + l) 

(3.2) 

On any given time slice, this global symmetry can be spon- 
taneously broken, and such breaking implies the absence of 
a confining color Coulomb potential. The color Coulomb po- 
tential can be extracted, at weak couplings, from the correlator 
of timelike links at a given time, i.e. |4] 



V c (/?) = -log 



\Ti[U (0,t)]U*(R,t)] 



(3.3) 



V C (R) converges to the instantaneous color Coulomb potential 
in the continuum limit. Asymptotically, this potential is also 
an upper bound on the static quark potential H12I1 



V(R) < V C (R) 



(3.4) 



so that a confining Coulomb potential is a necessary, but not a 
sufficient, condition for permanent confinement. 

The remnant symmetry breaking order parameter Q is ex- 
pressed in terms of the timelike link variables averaged, at a 
given time, over spatial volume (L 3 ) 



(3.5) 



If the remnant symmetry \2>.2\ is unbroken, then the modulus 
of U should vanish in the infinite volume limit, at any t . We 
therefore define the order parameter as 



1 L < 

l ' t=i 



{Tr[U(t)W(t) 



(3.6) 
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where L, is the lattice extension in the time direction. Q van- 
ishes in the large volume limit in the unbroken phase, and 
has a non-zero limit if the remnant symmetry is spontaneously 
broken. 

An exponential falloff in the timelike link correlator im- 
plies that the color Coulomb potential rises linearly with sep- 
aration. In contrast, if remnant symmetry is broken sponta- 
neously, then V C (R) — > constant as R — > °°. So the existence of 
an asymptotic Coulomb string tension G cou i > depends on 
the unbroken realization of remnant gauge symmetry, and this 
fact makes the Q order parameter a good candidate for isolat- 
ing the temporary confinement phase from the Higgs phase. 
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visible near y = 0.84. 3 Fig.Hshows our data for Q vs. y, again 
at (3 = 2.2, along with another observable s w to be discussed 
shortly. A sudden jump in Q is visible at the same value of 
y (within our resolution) that the jump in plaquette energy is 
observed. The non-zero value of Q below the transition is a 
finite-size effect; this quantity should vanish on an infinite lat- 
tice n. 



As p is reduced, the first-order transition disappears, as seen 
in the plot of plaquette energy vs. y at p = 1 .2, shown in Fig. 
|9] There is still, however, a transition in Q, as seen in Fig. 
ITUl The Q order parameter is not discontinuous in this case, 
instead, Q increases continuously away from zero (in the in- 
finite volume limit) upon crossing the transition line, at about 
y = 1.5. This behavior is reminiscent of magnetization in a 
spin system, in the neighborhood of a second order phase 
transition. In Fig. ^| we show data for both 8 4 and 16 4 lat- 
tices, to show the trend to Q — 0, at infinite volume, below the 
transition. We have determined the position of the remnant 
symmetry-breaking transition line at a range of couplings be- 
low p = 2.2; this is the lower line shown in Fig. II II 4 
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FIG. 7: Plaquette energy P vs. Higgs coupling y at (3 = 2.2; a weak 
first-order transition is seen. 
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FIG. 9: Plaquette energy P vs. Higgs coupling y at P = 1 .2; no tran- 
sition is evident. 



FIG. 8: Data for the remnant symmetry order parameter Q vs. y, and 
percolation order parameter s w vs. y, at (3 = 2.2. 



Further support for this view of Q as an order parame- 
ter comes from the fact that that the Q-transition line in the 
gauge-Higgs theory coincides with the thermodynamic first- 
order phase transition line, up to the terminating point of that 
line of transitions. Like the plaquette energy, the data suggests 
that Q is discontinuous along the transition line in the infinite 
volume limit. In Fig. we show the plaquette energy curve, 
as a function of y, at P = 2.2; a weak first-order transition is 



3 We cannot rule out the possibility that this is a very sharp crossover, rather 
than an actual first-order transition. 

4 The Q-transition line shown in Fig. II ll differs somewhat in location from 
the line we reported previously in Fig. 12 of ref. jQJ- The calculation of 
that figure suffered from an unfortunate program error; our current Fig. II II 
corrects and replaces it. We note that this program error did not affect the 
other results reported in ref. |4|]. 
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FIG. 10: At p = 1.2, the remnant symmetry and depercolation transi- 
tions occur at different values of y. Solid lines are presumed extrap- 
olations to infinite volume. 



B. Vortex Percolation 
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FIG. 11: Transition lines for the remnant symmetry breaking 
transition (lower line, order parameter Q) and the percolation- 
depercolation transition (upper line, order parameter ,s w ). The shaded 
region is a region of couplings where remnant symmetry is broken, 
but vortices still percolate. The two transitions lines appear to con- 
verge at the end-point of the line of first-order thermodynamic tran- 
sitions. 



We also have a second candidate for the role of order pa- 
rameter, distinguishing between the temporary confinement 
and Higgs phases. This is an operator which is sensitive to 
the vortex percolation-depercolation transition |3[, denoted 
s w , and defined as follows: Let f(p) be the fraction of the 
total number Np of P-plaquettes on the lattice, carried by the 
P-vortex containing the P-plaquette p. Then s w is the value of 
f(p) when averaged over all P-plaquettes. It can be thought of 
as the fraction of the total number of P-plaquettes on the lat- 
tice, contained in the "average" P-vortex. More precisely: let 
the index i = l,2,...,N v denote vortex number, and i(p) speci- 
fies the vortex containing the P-plaquette p. Also let denote 
the total number of P-vortices contained in vortex ;. Then 



N P p tf N P 



N * n 2 

L_L 
N 2 



(3.7) 



If all P-plaquettes belong to a single vortex, then s w = 1 . In 
the absence of percolation, the fraction of the total number 
of P-plaquettes carried by any one vortex vanishes in the in- 
finite volume limit. If a finite fraction of P-plaquettes is car- 
ried by a finite number of percolating vortices in the same 
limit, then s w > 0. The transition from s w > to s w = in the 
large volume limit identifies the percolation-to-depercolation 
transition. In ref. |3] we determined the line of depercolation 
transition in the gauge-Higgs model with variable Higgs mod- 
ulus; here we report the location (upper line in Fig. II 1> in the 
gauge-Higgs model J2.lt . for comparison with the remnant 
symmetry-breaking line (lower line in Fig. II 11 1. The calcula- 
tion of s w requires identifying, on each center-projected lattice 
configuration, the number of separate P-vortices and the area 
of each. Our algorithm for doing this is described in detail in 
the Appendix. 

The vortex depercolation transition s w — » 0, like the Q- 
transition, coincides with the line of thermodynamic, first- 
order transitions up to the endpoint of that line. The data for 



s w vs. y at p = 2.2 is displayed, together with the Q data, in 
Fig-IU However, beyond the first-order transition line, the vor- 
tex depercolation and remnant symmetry-breaking transitions 
no longer coincide, as is evident from our data for Q and s w at 
P = 1.2, displayed in Fig.^5| At this coupling, the Q transi- 
tion occurs at about y = 1.5, while s w goes to zero at y = 2. 1 . 

The transition lines for both remnant-symmetry breaking 
and depercolation can be compared in Fig. ^2 over a range 
of couplings p, and it is evident that these transition lines do 
not coincide, contrary to what was assumed implicitly in ref. 
0| . There is a region between the two transition lines where 
vortices percolate, but G cou i — 0. This is consistent with the 
notion that vortex percolation is a necessary, but not sufficient, 
condition for having a confining Coulomb potential, which is 
itself a necessary, but not sufficient, condition for permanent 
confinement. 

The fact that the vortex and remnant-symmetry transition 
lines do not coincide tends to support the most straightfor- 
ward interpretation of the Fradkin-Shenker theorem, namely, 
that there is no unambiguous distinction between the tempo- 
rary confinement and Higgs phases. In either region, the large- 
scale gauge-field fluctuations responsible for disordering Wil- 
son loops are suppressed, and the gauge field due to an exter- 
nal static source falls off exponentially with distance from the 
source. In this sense the regions are very much alike in the 
far infrared, and are characterized by charge screening, rather 
than confining forces. 



IV. CONCLUSIONS 

There are two conclusions. First, the vortex mechanism for 
producing a linear potential can work even when the gauge ac- 
tion does not possess a global center symmetry, and the static 
potential is flat at large distance scales. Thus global center 
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symmetry is not necessarily essential to the vortex mecha- 
nism. The same distribution of vortices which produces a 
linear potential over a finite interval, in temporary confine- 
ment theories, can also avoid producing a linear potential at 
asymptotic distances, as seen in the Polyakov line correlator 
on the center projected lattice. Of course, the distribution of 
percolating P-vortices responsible for permanent confinement 
at y = 0, and that responsible for temporary confinement at 
y > 0, must differ qualitatively in some way at large scales. In 
the latter case, we would expect that vortex piercings of the 
minimal area of a very large Wilson loop would tend to come 
in pairs, whose effect on the large loop would cancel. Whether 
this effect is due to vortices having a branched polymer struc- 
ture at large scales, or is due to some other distribution, is left 
for future investigation. 

The second conclusion concerns the question of whether it 
is possible, in a theory without a local order parameter for 
confinement, to nonetheless distinguish between a "confined" 
phase and a Higgs phase via some non-local order parame- 
ter. We have investigated two reasonable candidates: (i) the Q 
observable which tests for spontaneous breaking of remnant 
gauge symmetry in Coulomb gauge, corresponding to the loss 
of a confining color Coulomb potential, and (ii) the s w ob- 
servable, which tests for P-vortex percolation. Both of these 
observables are closely related to confinement in pure gauge 
theories; unbroken remnant symmetry is a necessary condi- 
tion for permanent confinement, and vortex removal removes 
the confining force. Moreover, both observables have a transi- 
tion in the gauge-Higgs phase diagram which agrees with the 
first-order transition line, up to the endpoint of that line. Be- 
yond the first-order transition line, however, we find that the 
remnant-symmetry breaking and vortex depercolation lines do 
not coincide, which means that the separation of the gauge- 
Higgs phase diagram into a "confinement" phase and a Higgs 
phase is ambiguous. The choice of a particular non-local ob- 
servable to be an order parameter for confinement is not very 
compelling, if the only non-analytic behavior seen at the tran- 
sition is in that particular observable. The fact is that through- 
out the phase diagram, the gauge-Higgs model at large scales 
is best described as a color-screening phase. In this model 
there are no large-scale gauge field fluctuations, characteristic 
of confinement, which disorder Wilson loops, and the color 
field due to a static source is screened (as in an electrically 
charged plasma, or in an electric superconductor), rather than 
collimated into a flux tube. This observation, together with 
our numerical result, tends to support the most straightforward 
reading of the Fradkin-Shenker theorem: There is no essential 
distinction, in a gauge-Higgs model, between the temporary 
confinement phase and the Higgs phase. 
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APPENDIX 

In this appendix we describe our procedure for identifying 
individual P-vortex surfaces. The basic idea is that two P- 
vortex plaquettes on the dual lattice which share a common 
link must belong to the same P-vortex. An ambiguity arises 
when four or six plaquettes share a link. This could be a self- 
intersection of a single P-vortex, or an intersection of two or 
more separate P-vortices. We simply ignore these ambiguous 
links; they are not used to identify different plaquettes as be- 
longing to the same vortex surface. The algorithm goes as 
follows: 

I. Gauge fix the SU(2) lattice configuration to maximal center 

gauge, and center project. The center-projected plaque- 
ttes have values z M v(x) = ±1, where x is the lattice site, 
and /j, v specifies the plane of the plaquette. 

II. Map the above plaquette variables onto variables onto pla- 

quette variables of the dual lattice 

£p{x+-jl+v)=zuv{x) (A.l) 

Although by convention the coordinates of points on the 
dual lattice are half-integer, we have added a constant 
vector to all dual lattice sites in order have 

integer coordinates on the dual lattice also. 

III. Count the number of negative plaquettes on the dual lat- 

tice, and assign to each of these a number « a p(x) from 1 
to A, where A is the total number of negative plaquettes. 

IV. Initialize n/ = 0. Loop through all of the links of 

the dual lattice. For each link shared by two and 
only two negative plaquettes, increment «/, and store 
the plaquette numbers of the two negative plaquettes 
in p(ni, 1), p(ni, 2). We will refer to such links as 
"surface-pair" links. Upon completion of the loop over 
links, set Ni equal to the final value of «/; this is the total 
number of surface-pair links. 

V. Initialize the Vortex Number of each negative plaquette, 

V(n) —0, n — l,...,A, and set n v = 0. Now loop through 
surface -pair links, n/ = 1, ...,Ni. At each such link de- 
note p\ = p(ni,l),p2 = /?(«/, 2), and then perform the 
following operation on the vortex numbers: 

1. if V(pi) = V(p2) = 0, increment n v — » n v + 1, and 
setV(p 1 )=V(p 2 )=n v . 

2. ifV(pi)^0, V(p 2 )=0, setV(p2)=V(pi). 

3. ifV{p 2 )^0, V(pi)=0, setV(pi)=V(p2). 

4. if both V(pi), V(p2) are non-zero, and V(p\) < 

V{ P 2),SCtV(p2) = V{ Pl ). 

5. if both V(pi), V(p 2 ) are non-zero, and V(p2) < 
V(pi),setV(pi) = V(p2). 
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6. if both V(pi), V(p2) are non-zero, and V(p\) = 
V(p2), do nothing. 

At the end of looping through the surface-pair links, set 
N v = n v . 

VI. Repeat step V, except for the initializations and the setting 

of N v . Continue iterating through the surface-pair links 
until convergence is reached; i.e. there is no further 
modification of the {V («)}. 

VII. Initialize vortex areas a„ = 0, n = 1,...,N V . Loop 
through the negative plaquette number n = 1,...,A. At 
each plaquette, increment 

a V(n) = a V{n) + 1 (A.2) 
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